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The heat generation by an electric current flowing through a quantum dot with the dot containing both 
electron-electron interaction and electron-phonon interaction, is studied. Using the non-equilibrium Keldysh 
Green's function method, the current-induced heat generation is obtained. We find that for a small phonon fre- 
quency, heat generation is proportional to the current. However, for a large phonon frequency, heat generation is 
in general qualitatively different from the current. It is non-monotonic with current and many unique and inter- 
esting behaviors emerge. The heat generation could be very large in the Coulomb blockade region, in which the 
current is very small due to the Coulomb blockade effect. On the other hand, in the resonant tunneling region, 
the heat generation is very small despite a large current, an ideal condition for device operation. In the curve of 
heat generation versus the bias, a negative differential of the heat generation is exhibited, although the current 
is always monotonously increasing with the bias voltage. 
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I. INTRODUCTION 

With the continued miniaturization of integrated circuits, 
quantum phenomena and new technological problems are 
emerging quickly. An essential technological roadblock of 
this miniaturization is the "power problem". With millions 
of transistors assembled on a chip area no larger than a few 
square centimeters, the chip-level power densities may reach 
to the order of lOOW/cm'^. The large amount of accumu- 
lated heat caused by the current makes the chip's temperature 
rise to a level so high that the chip may not function properly. 
There are two ways to reduce the chip's temperature: one is 
to remove the heat as quickly as possible, while the other is to 
suppress the heat generation. It is unfortunate that the novel 
and complex device geometries normally make heat removal 
more difficult and most of the new materials appearing in de- 
vice processing have lower thermal conductivities than bulk 
silicon.^ Thus, it is essential to uncover the laws of heat gen- 
eration induced by an electric current in nanoscale devices and 
to investigate how the heat generation may be reduced. 

The idea that electric current does work on conductors and 
causes heating has been known since the mid-nineteenth cen- 
tury. It has been well studied in macroscopic systems and 
is known as the effect of joule heating; however, many open 
questions remain concerning heating in nanoscale systems. 
For a net resistive nanoscale system, the total heat dissipa- 
tion is still equal to the input power IV, with / being the cur- 
rent / and V the bias voltage. In standard transport theory, 
such as the Landauer-Buttiker formula or the non-equilibrium 
Keldysh Green function (NEKGF) method, it is normally as- 
sumed that the electron is elastically scattered in the central re- 
gion (called the scattering region) and that all heat dissipation 
occurs in the reservoirs.^ This assumption does not provide a 
full answer to the dissipation problem. One can see this from 
the following two aspects: (1) If the central region contains 
the phonon degrees of freedom and the electron-phonon (e-p) 
interaction, the inelastic scattering processes and heat dissipa- 
tion can occur in the central region. In this case, the simple 



assumption of all dissipation occurring in the reservoirs is in- 
valid. In fact, in some recent experiments, it is found that the 
e-p interaction is indeed strong in some nanoscale molecular 
devices, such as. Ceo, nanotube, etc.^ ^^'^ In these devices, the 
phonon-assisted tunnelling sub-steps or sub-peaks have been 
observed in the I-V curves or in the differential conductance 
versus the gate voltage. ^ "^'^'^'^ (2) Furthermore, if one wishes 
to investigate where and how the heat dissipation occurs while 
a current passes a nano-device even assuming the device con- 
tains no e-p interaction, the e-p interactions in the leads that 
near the nano-device must also be considered. Then this part 
(i.e. the dissipation region) can not be considered a part of the 
reservoirs since the reservoirs are assumed to be without any 
interactions in the standard transport theory. 

Recently, some attention has been paid to the problem of 
current-induced heat generation in nanoscale devices, with a 
few preliminary results. ^^ ^^^^i^^i^^i^^ Huang et al have exper- 
imentally observed the current-induced local heating effects 
in single molecules by measuring the force required to break 
the molecule-electrode bonds Very recently, Oron-Carl and 
Krupke have successfully determined the hot phonon gener- 
ation from a current bias by using the ratio of anti- Stokes 
to Stokes lines Theoretically, Lazzeri et al have stud- 
ied the heat generation in nanotube via the first principles 
calculations.^^ All previous results demonstrate that the heat 
generation in nanoscale devices is important. In a recent 
work, we studied the heat generation when an electric current 
passes through a mesoscopic devices. By using the NEKGF 
method, a general formula for the current-induced heat gen- 
eration was derived. This formula can be applied to both the 
linear and nonlinear transport region, as well the mesoscopic 
device with various interactions. 

In this paper, we study the current-induced heat generation 
while an electric current flows through a QD that is in the 
Coulomb blockade regime, focusing on the heat transfer to 
the phonon part since the phonon creation is the main source 
of the heat generation in QD. Considering that the device con- 
sists of a QD coupled to two leads and the QD contains the 
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FIG. 1: (Color online) The heat generation Q(Ag) and current / vs. 
the bias V for the different temperature T with the level = (a) 
and 1 (b), and the other parameter is cjg = 3, [/ = 20, Ag = O.Gcjq. 
The inset in (b) is schematic diagram for the tunneling process in 
which a phonon is absorbed. 

electron-electron (e-e) Coulomb interaction and e-p interac- 
tion, by using the NEKGF method and the general formula 
of heat generation, both the heat generation and the current 
are obtained. The numerical results exhibit that the heat gen- 
eration in the Coulomb blockade regime has many interest- 
ing and unique behaviors. For example: (1) In the resonant 
tunneling region with the intra-dot level within the bias win- 
dow, the current is very large but the current-induced heat 
generation is very small. This is an ideal region for the op- 
eration of the nano-device with a large current and low heat 
generation. On the other hand, in the Coulomb blockade re- 
gion, the current is very small but the heat generation is very 
large. In particular, when the e-e interaction strength is equal 
to the phonon frequency, a resonant phonon emitting process 
occurs and a large amount of current-induced heat genera- 
tion emerges, even though the current is quite small due to 
the Coulomb blockade effect. (2) In some parameter regions, 
when the bias rises, the current increases but the heat genera- 
tion decreases. In other words, the negative differential of the 
heat generation emerges. 



II. MODEL AND FORMULATION 

The system of the lead-QD-lead can be described by the 
following Hamiltonian: 

a,k,a a,k,a 
a 

where fia- = d\dcr, and a = L^R represent the left and right 
leads, cj^^^ and dj. create an electron with spin a in the a lead 
and QD, respectively. Analogously, aj is the phonon creation 
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FIG. 2: (Color online) The of heat generation Q(A^) vs. the phonon 
frequency cuq for different e-e interaction strength U with the param- 
eters: y = 4, T = 1, Ag = O.Gcjg, and = 

operator and ujq is the phonon frequency. The QD has a sin- 
gle energy level Sd with the spin a =T, j. The intra-dot e-e 
Coulomb interaction and e-p interaction are considered, with 
the interaction strength U and A^. The second term in Eq. ([T]) 
describes the tunneling coupling between the QD and the two 
leads and t^k is the hopping matrix element. 

Before to solve the heat generation in the QD we first 
make a canonical transform with the unitary operator^^ U = 
exp{— ^cr{^q/^q){dl — dq)ha}. After this transformation, 
the Hamiltonian becomes: 

H = ^edd^da + Un^fll + ^ Sakclka^aka 
c ak,a 

+ ^ {takcl^daX + H.c.) + ujqa\aq, (2) 

where £d = €d - )^l/^q, U = U — 2\yuOq, and 
X = exp{ — (Ag/cjg)(aJ — dq)}. We make an ap- 
proximation to replace the operator X by its mean value 
(X) = e^v{-{\/^q?{Nq + 1/2)}P where Nq = 
l/[exp{uq/kBT) — 1] refers to the phonon number. This ap- 
proximation was used in previous studies,^' and it is known 
to be valid when tak ^q- After taking this approximation 
the phonon is decoupled with the electron. Then the formulate 
of the current-induced heat generation Q per unit time is'i^ 

-2N,[G>,icu)GU^) + Gl,(co)G>,im- (3) 

where O = co — coq, and (5^'^'^'^(cc;) are the Fourier trans- 
forms of G^'^'^'^(t), which are the standard QD's single- 
electron Green's functions of the Hamiltonian (O and are 
defined as G^,(t) = -ie(t)({d(i), dt(0)}), G-,(t) = 
iei-t){{dit),d^iO)}),G<At) = i{S{Q)d{t)),G>,{t) = 
—i{d{t)d) {{))). Here heat generation Q is the dissipation in 
the QD region in which the energy is transferred to the local 
phonons. Since the present device is a net resistive system, the 
total dissipation is equal to the input power IV, and remaining 
dissipation Qres occurs in the reservoirs. Moreover the heat 
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III. NUMERICAL RESULT 
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FIG. 3: (Color online) Q{Xq) and / vs. the level id for the different 
U with the parameters a;^ = 8, F = 4, T = 1, and — 0.6oOq. 
The insets in (b) are the schematic plots for two phonon-emitting 
processes. 

generation actually is the rate of heat generation. Hereafter we 
don't distinguish them because that they have the completely 
same characteristic. Because of the e-p interaction has been 
decoupled in the Hamiltonian Q, the Green's function G^^ 
can be obtained by the equation-of-motion technique i^ii 



U{1 - (na))] 



(4) 



where oje = oj — id, (n^) is the electron occupation num- 
ber in the intra-dot state <?, a =[ while a =] and a =] 
while a =lt = (fi + ffl)/2, f« = r«(X)2, and 
= 27r|tttfcp(5(w — Cak) is the linewidth functions 
which assume to be independent of the energy tu?-^'' Af- 
ter solving G^,(c^), = G<,{u;) = 

G^{u;)t<{u;)G-{u;l and G^ = G<, + G^, - G^,, 
where t<{u;) = ilTLfUu;) + tRfniu;)] and U{u;) = 
l/{exp[{uj — /jja)/kBT] + 1} is the Fermi distribution func- 
tion in the a lead. Substituting these Green's functions into 
the Eq.©, the heat generation Q can be obtained straight- 
forwardly. Finally, the electron occupation numbers (ricr) in 
Eq.© need to be self-consistently calculated with the equa- 
tion (ria) = i j {duj /2it)G^^{<jj). In addition, by using the 
same method as in Ref.^, the current can also be calculated. 



In the numerical investigation, we consider the symmetric 
barriers^^ with V l = T ^ = T, and set F = 1 as the unit of 
energy and /ii? = as the energy zero point, then the bias 
eV = fiL — fJ^R = Ml- Fig. la shows the heat generation Q 
and current / as functions of the bias V for the different tem- 
perature T at the QD level id = 0. We can see that at first the 
heat generation Q is almost zero while the bias V < ujq since 
the bias can not offer enough energy to emit in this range. 
With the bias passing coq, both of Q and / rapidly increase 
and reach their respective plateaus. Q and / remain at plateaus 
until the bias reaches U, then they rapidly increase and reach 
new plateaus. The reason behind the two steps is that the QD's 
level is split into two sub-levels at id and id^U with the exis- 
tence of the e-e interaction U. The temperature has little effect 
on the heat generation. Raising temperature makes the curves 
of both Q-V and I-V more smooth, but the qualitative behav- 
iors remain. Here we emphasize that the heat generation Q is 
not directly proportional to the current /, and they obviously 
has the following two differences: (i) Many small sub-steps 
emerge in the current curves because of the phonon-assisted 
tunneling processes, but no phonon-assisted sub-steps in the 
heat generation curves, (ii) The rapid jumps in the current 
curve happen at about V = id and id + but the rapid jumps 
in the heat generation curve happen at about V = id^ojq and 
id^U ^ ujq, a delay of coq (see Fig. la). 

In Fig. la, the level id is fixed at zero. Next, we study 
the case of id 7^ 0. Fig. lb describes the heat generation Q 
with the same parameters as in Fig. la except q = 1. The 
results clearly exhibit that a negative differential heat gener- 
ation dQ I dV emerges when V is between U -\- id — ujq and 
U ^ id. In this region, with increasing V, I increases, but 
Q decreases. The reason is that the electron can absorb as 
well as emit a phonon while it tunnels through the QD. When 
the bias V is between U -\- id — ooq and U -\- id, the tunnel- 
ing process as shown in the inset of Fig. lb occurs, in which 
the electron absorbs a phonon while it passes through the QD. 
Thus, a negative differential heat generation emerges. 

Fig. 2 investigates the dependence of the heat generation Q 
on the phonon frequency ujq. When ujq approaches to zero, so 
does the heat generation Q. With coq increasing from zero, Q 
quickly increases. The curve of Q-cOq exhibits two peaks. One 
peak is roughly atcOq = T, very close to zero. The other peak 
is at cc^g = U. The first peak is caused by competition of two 
effects when raising coq'. (1) the electron losses more energy 
when it emits a phonon, at the same time (2) it is more and 
more difficult to emit a phonon. The second peak is largely 
associated with the e-e interaction strength U. Its location is 
fixed at coq = U and its height strongly depends on U. The 
smaller U is, the higher the peak is. In fact, this peak origins 
from a resonant phonon emitting process as shown in the right 
inset of Fig. 3b, which will be discussed later in more detail. 

Fig. 3 shows Q and / versus the QD's level id for different 
values of the e-e interaction U. I has two main resonant peaks 
at id = p and id-\-U = with /i = (/iL+Mi?)/2. In addition, 
some phonon-assisted sub-peaks also emerge at the two side 



4 



0.40 0.80 1.2 1.6 20 0.40 0.60 0.80 




-9 -6 -3 3 6 -9 -6 -3 3 6 

FIG. 4: Q(Xq) (a) and / (b) vs. the level id and phonon frequency 
ouq with = 6, y = 4, T = 1, and Ag = O.Guq 

of the main peaks, and the distance to its main peak is ±nujq 
(n = 1,2,3,...). These well-known results on current have 
been experimentally observed^ Following, we focus on the 
heat generation Q. The curves of Q versus the level id have 
three peaks, which are marked 'A', 'B', and 'C. The peak 'A' 
('B') stands at the valley between the right (left) main peak 
and its right-side (left- side) sub-peak of the current. The peak 
'C places at id -\- U /2 = j2, i.e. at the valley (the Coulomb 
blockade region) between the two main peaks of the current. 
Here we emphasize that none of the peaks of Q aligns with 
that of /, in fact all peaks of Q appear at the valleys of /. 
This means that the heat generation is small for large current, 
and vice verse. Here we also emphasize another important 
feature: with the e-e interaction U approaching to the phonon 
frequency ujq, the peak 'C is greatly enhanced. AtU = ujq, 
the peak 'C is very large. For example, the peak 'C in Fig. 3b 
with U = ujq = 8 is hundredfold higher than the peak 'C in 
Fig. 3a with U = 10. Due to the largeness of peak 'C, the 
peaks 'A and 'B' are not visible in Fig. 3b, although they are 
there. As far as for the current, there is no such high peak at 
all. 

Let us explain why the heat generation Q is large at the 
valleys of the current. This is because the phonon-emitting 
processes are stronger in the current valley region. For ex- 
ample, to consider the phonon-emitting tunneling process in 
the left inset of Fig. 3b, in which an electron from the left lead 
tunnels into the QD at the level id, emits a phonon and jumps 
to the sub-level id — ujq due to the e-p interaction, then tun- 
nels to the right lead. The strength of this process is direct 
ratio of the occupation probability fhi^d) of the left lead at 
the energy co = id, and the empty probability I — fnicd — ujq) 
of the right lead at the energy co = id — ujq. So this process 
is maximum at id — ujq/2 = ft, and leads to the peak 'A in 
the curve of Q-id- The origin of the peak 'B' is same as for 



the peak 'A. The peak 'C origins from the tunneling process 
as shown in the right inset of Fig. 3b, which is similar to that 
in the left inset of Fig. 3b. But the two levels id and id ^ U 
are real electronic levels now. Therefore the resonant phonon 
emitting occurs when the difference 0) of two levels is equal 
to the phonon frequency ujq, giving rise to a huge peak in heat 
generation. 

Finally, in Fig. 4 we present the contour plot of the heat gen- 
eration Q and current / versus the QD level id and phonon 
frequency ujq. As Uq approaching to zero, Q also approaches 
to zero regardless of the other parameters, however, / can be 
large if the resonant tunneling occurs. At small but non-zero 
cjg, e.g. ujq < r, the heat generation Q is approximatively 
proportional to the current /. In this case the quantum re- 
sult of the heat generation resembles the classical one. On the 
other hand, at large ujq (e.g. ujq > F), the contour of the heat 
generation is totally different to the contour of the current. For 
example, for uoq = 17,3. great heat generation appears, but we 
can not find any large change in the current plot. 



IV. CONCLUSION 

In summary, we investigated the heat generation in a quan- 
tum dot (QD) under a bias and in the Coulomb blockade 
regime. The results show that for the small phonon frequency, 
the heat generation is proportional to the current. On the other 
hand, for the large phonon frequency, the heat generation is 
quite different from the current and exhibits many unique and 
interesting characteristics. For instance, a negative differen- 
tial heat generation emerges in some parameter regions even 
though the differential conductance is always positive. In par- 
ticular, the heat generation is large in the Coulomb blockade 
region, while it is small in the resonant tunneling region. The 
later behavior is advantageous for the QD device because the 
QD device normally works in the resonant tunneling region 
and low heat generation is ideal for device operation. 
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